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Abstract 

In this paper we extend the conditions on quasi-thermal-incompressible materials 
presented in [1] so that they satisfy all the principles of thermodynamics, including 
the stability condition associated with the concavity of the chemical potential. We 
analyze the approximations under which a quasi-thermal-incompressible medium 
can be considered as incompressible. We find that the pressure cannot exceed a 
very large critical value and that the compressibility factor must be greater than a 
lower limit that is very small. The analysis is first done for the case of fluids and 
then extended to the case of thermoelastic solids. 
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1 Introduction 



It is well known that fully incompressible materials do not exist in nature. 
However, it is important to have a mathematical model for an incompressible 
medium, as an idealization of media that exhibit extreme resistance to volume 
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change. In the context of isothermal mechanics, an ideal incompressible ma- 
terial is a medium that can only be deformed without any change in volume. 
Extensive literature has been devoted to qualitative analysis and numerical 
methods for constructing solutions of incompressible fluids as limits of com- 
pressible ones as the Mach number tends to zero (see for example [2|3] ). 
However, when the process is not isothermal, the notion of incompressibility 
is not well defined. Several possibilities arise. 

For compressible fluids, the pressure is a constitutive function while for incom- 
pressible fluids the pressure is only a Lagrange multiplier associated with the 
constraint of incompressibility. Therefore, to compare the solutions of com- 
pressible and incompressible media, it is convenient to choose the pressure p 
(instead of the density p) and the temperature T as thermodynamic variables 
(see, e.g. [i|5] ). the other quantities, such as specific volume V = 1/p and 
internal energy e, being determined by constitutive equations in the form: 

V = V(p,T), e = e(p,T). 

Two parameters are important for a fluid: the thermal expansion coefficient a 
and the compressibility factor (3 defined by 

Vt r V p m 

a = v> /3 = _ y (1) 

where the subscripts T and p indicate partial derivatives with respect to vari- 
ables T and p. 

Experiments confirm that for fluids considered as incompressible the volume 
changes little with the temperature and remains practically unchanged with 
the pressure. For this reason, many authors consider as incompressible a ma- 
terial for which the specific volume does not vary with the pressure but varies 
only with the temperature (i.e., V = V(T)). 



The first model of incompressibility was proposed by Miiller [6]: Here, all 
the constitutive equations of an incompressible fluid do not depend on the 
pressure. We think that Midler's motivation stems from his using of variables 
p and T: in experiments, the density of incompressible materials depends on 
the temperature and it is reasonable to assume that this is the case for all 
constitutive functions as for instance for the internal energy. 
Nevertheless, Miiller proves that the only function V = V(T) compatible 
with the entropy principle is a constant. As previously indicated, this result 
obviously disagrees with experimental or theoretical results as in the so-called 
Boussinesq approximation (see, e.g. [TJ, [8]). We called this contradiction the 
Miiller paradox pp. 

A second, less restrictive, model usually employed in the literature requires 
that the only constitutive function independent of the pressure is the specific 
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volume (see, for example Rajagopal et al [US]) and as a result the Gibbs 
equation is satisfied. In a recent paper [I], we named such a material a quasi- 
thermal-incompressible medium and proved that for a pressure smaller than a 
critical value, the ideal medium of Miiller can be recovered as a limit case. 
Nonetheless, a weakness of V = V(T) as definition for incompressibility was 
first noted by Manacorda [9] who showed that instabilities occur in wave prop- 
agation. The instabilities are due to the non-concavity of the chemical poten- 
tial and to the sound velocity c becoming imaginary; the mathematical system 
of Euler equations is of elliptic type. Let us also note that in the solid case, 
Dunwoody and Odgen [12J showed that conditions for infinitesimal stability 
are unattainable if the trace of the strain tensor is dependent on tempera- 
ture. Starting from Manacorda's observation, other authors such as Scott et 
al [TOfTI] proposed an alternative definition of incompressibility: instead of 
V = V(T), they assumed V = V(S), where S is the specific entropy; then, in 
the ideal incompressible case, the Mach number becomes unbounded and the 
sound velocity becomes infinite. From an experimental point of view, this as- 
sumption is unrealistic because the entropy is not an observable and no direct 
evaluation of V = V(S) is possible. Moreover, this ideal limit case implies a 
parabolic structure for the Euler equations. 

Due to the above considerations, the aim of our paper is to propose a more 
realistic model of incompressible medium as limit case of a quasi-thermal- 
incompressible material based on the following requirements: 

(i) The model must respect all the principles of thermodynamics: both the 
Gibbs equation and the thermodynamic stability corresponding to the con- 
cavity of the chemical potential must be satisfied. 

(ii) The model must fit with experiments when the compressibility factor is 
not zero but it is very small. 

It is noteworthy that we obtain numerical results perfectly fitting with those 
obtained in paper pQ . The results are also extended to hyperelastic materials. 



2 Thermodynamic restrictions 

We consider the thermodynamic conditions verified by compressible fluids 
when the specific volume is governed by a constitutive function written in 
the form: 

V = V{p,T). (2) 
The entropy principle and the thermodynamic stability must be satisfied. 

a) The entropy principle 

In local equilibrium, the entropy principle requires the validity of the Gibbs 
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equation: 

TdS = de + pdV. (3) 

In addition for Navier-Stokes-Fourier fluids, the heat conductivity and the 
viscosity coefficients must be non negative although they are null for Euler 
fluids. The choice of independent variables p and T induces the chemical 
potential pu clS db natural thermodynamic potential: 

l i = e+pV -TS 

and Eq. ([3]) is equivalent to 

dp, = Vdp — SdT. 
Other thermodynamic variables derive from the chemical potential: 

V = f i p , S = -fi T (4) 
and the specific internal energy is 

e = [i - p p p - T fpr, (5) 
where for any function / = f(p,T), we denote 

MIL /r= (S) P - 

The thermal equation of state (j2J) is determined by experiments while the 
chemical potential and the entropy density can be deduced as follows: 

H= Jv(p,T)dp + jl(T), S=-Jv T (p,T)dp-jl'(T), (6) 

where Jl (T) is a function only depending on T. From Eq. (jSJ) we get 

e = e (T) + Jvdp-pV -T J V T dp, (7) 



with 



e(T) = Jt{T)-Tp, '(T) or equivalently Ji(T) = -T J dT. (8) 
We summarize this statement: 

Statement 1 - For any constitutive functions V = V(p,T) and e = e(T), 
the entropy principle is satisfied if the chemical potential, entropy density and 
internal energy are given by Eq. (jHJ)i, Eq. ([6]) 2 and (JTJ), together with Eq. flEJta- 

b) Thermodynamic stability 
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The specific heat C P is defined as the partial derivative of the specific enthalpy 
h = e + p V with respect to T at constant pressure p. Consequently, Eq. (J7J) 
yields 

C P = e'(T) — T J V TT dp . (9) 

Thermodynamic stability requires that the chemical potential be a concave 
function of p and T: 

Vpp = v p < 0, U TT = ~ < o) (10) 

and 

C V TV 2 

i = fi TT fi pp - fi 2 TP = — -v^>o <^> v v < (ii) 

By using Eq. ([T]), inequality ( II ip can be written in terms of the thermal 
expansion coefficient a and the compressibility factor (3: 

a 2 TV 

> Per, Per = -pr- > 0. (12) 

Up 

Statement 2 - Thermodynamic stability requires that the state functions V 
V [p. T) and e = e(T) satisfy the inequalities: 

TV 2 

Vp < — pr^, Cp > 0. 
Up 

Consequently, there exists a lower bound limit (3 cr of ft such that if (3 > Per, 
then the material is stable. 

The adiabatic sound velocity c is: 
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From Eq. @ we obtain 

dV = fipT dT + fi pp dp, dS = ~Htt dT — fi Tp dp. (13) 

When dS = 0, Eq. (JIB"])? substituted in Eq. f TT5]) 1 yields p as a function of V 
and we get: 

c 2 = (14) 

Therefore, when the chemical potential is a concave function of p and T, 
we automatically get c 2 > and the differential system for Euler fluids is 
hyperbolic. Taking account of Eqs. ©1, ©2, ©, (HU) and flTJ), Eq. JUj) 
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yields the sound velocity in the form: 




We observe that when V = V(T) the model of incompressibility corresponds 
to j3 = and for Euler fluids the differential system is elliptic; when V = 
V(S), the model of incompressibility corresponds to = /3 cr and the system 
is parabolic. In our case /3 > (3 C r, the differential system is hyperbolic and the 
fluid is stable. 



3 Quasi-Thermal-Incompressible Materials 

For the so-called incompressible fluids, the volume changes little with the tem- 
perature and changes very little with the pressure. Experiments confirm this 
assumption. 

In the neighborhood of a reference state (p ,T ,V ), we choose a small dimen- 
sionless parameter 8 (S «C 1) such that: 

5 = a T , (15) 

and moreover, we assume that (3 is of order 5 2 : 

(3 oPo =0(5 2 ), (16) 

where a and /3 are the thermal expansion coefficient and the compressibility 
factor at the reference state. 

Definition 1 - A compressible fluid satisfying the thermodynamic conditions 
of Section 2 is called an Extended- Quasi-Thermal-Incompressible fluid (EQTI) 
if there exist V(T) and e(T) such that 

V{p,T) = V{T)+0{5 2 ) with V'(T) = 0{5) and e{p,T) = e{T)+0{5 2 ). (17) 

This means that an EQTI material is a stable compressible fluid that approx- 
imates an incompressible fluid to order 5 2 in the sense of Muller's definition. 
Conditions f|T5l) - f|T6l) together with Eq. f|T7|) i yield the representation of V(p, T): 

V(p, T) = V + 8 W(T) - 8 2 U(p, T), (18) 

where W(T) and U(p,T) are two constitutive functions chosen in agreement 
with conditions in Section 2. From Eqs. (jn])i and (ITS!) we deduce: 

H = ii(T)+pV + 6p W(T) + 8 2 fi(p, T), fi(p, T) — — J U(p, T) dp. 
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Since we can incorporate its limit value in \l (T) without loss of generality, we 
can chose the function p, so that lim p „ >0 p,(p, T) = 0. From Eq. (J7|) we obtain: 

e(p,T) = e(T)-5TW'(T)p + 0(5 2 ). (19) 

Due to Eq. ffT9l . in order to satisfy Eqs. f lT7|) 9 we require that the pressure 
cannot exceed a critical value: 

p« Pcr (T) with p cr {T) = - 5 T e ^ T) . (20) 

We observe that the critical value is large and of order 8 . In incompressible 
materials, the pressure cannot exceed a critical value depending on tempera- 
ture. Inequality (|20|) implies 

e(p,T) = e(T) + 0(5 2 ). (21) 

Moreover, from Eqs. ([1]) and (fT5l) we obtain 

a = d —(7 — + 0{°)> w ( T o) = ^t; p = ° — t? — ■ 

^0 ' 



From Eqs. fl9J) and (jT2j) we get 



. 2 rW(r) 



/3 cr = 5 2 — — ^, C P ~e'(T), 
K Up 



and the thermodynamic stability is ensured when 

TW' 2 (T) 

U p (p, T) > { \ C P ~ e'(T) > 0. (22) 

Finally, from Eq. ( TT4l) we deduce that the dominant part in sound velocity is 
of order <5 _1 : 

8\U P C»-TW' 2 ■ [ ] 

We can conclude with 

Statement 3 - An EQTI fluid given by constitutive functions f fl8|) and (12 ip . 

satisfying inequality fl22|) zs a g'oorf approximation of an incompressible fluid as 

V and e differ to order 5 2 from functions depending only on T, provided the 
pressure is smaller than a critical pressure p cr given by Eq. fl2TJ|) . The sound 
velocity given by Eq. f[2Uj) is real. 

Remark: We notice that in the limit case of isothermal processes for which 
W(T) = and U = U(p), we have (3 cr = and inequality ( )22|) i and the sound 
velocity become U'(p) > and c = Vo/^Sy/TF), respectively. In this case again, 
the EQTI requires that V is not constant but is function of p in the form 

V = V -5 2 U(p). 
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3. 1 Linear dependence of V with respect to T and p 



The most significant case is the linear expansion of V near (T ,p ): 

V = V {1 +a(T-T ) - P(p-p )} with e = C P T. (24) 



In this case, the scalars a, (3 and C P are positive constants. Expression (|24|) is 
a particular case of Eq. ( jl8j) under the identifications: 

aT =5, W(T) = ^(T-T ), U(p,T) = £^ (p - p ). (25) 



Then, the fluid is EQTI if 



aT = 5<l and /3 > (3 cr with /3 cr = 5 2 -^-, (26) 



together with 

Relation ( 1231) and relation (|25l) yield the adiabatic sound velocity: 



1 C T 

p<.p cr , with p cr = 7 — j— 9 -. (27) 

5 V 



/3 C P - a 2 V T C 2 C P 



Consequently, a fluid can be considered as incompressible if the pressure is 
smaller than a critical pressure which is of order a -1 . This is similar to con- 
clusions in paper pQ. In addition, to ensure the convexity of the chemical 
potential, the compressibility factor must be very small but not identically 
null: P must be greater than a critical value /3 cr which is of order a 2 . From 
Eq. (|28|) . the value of the sound velocity allows to calculate the value of j3. 
We finally observe that, in the Boussinesq approximation, (3 ^ induces an 
additional term depending on p in the density expansion near (T ,p ): 

p = p {l -a(T-T )+ (3{P~Po)}- 

Therefore, these considerations should be useful to revisit and to justify more 
rigorously the Boussinesq approximation. We point out that in this case we 
obtain constitutive equations that may be implicit in the direction pointed 
out in 



3.2 Case of liquid water 

Numerical values for liquid water are obtained in [14j : 

p = 10 5 Pascal; V = 10" 3 m 3 /kg; T = 293 °K; C P = 4.2 x 10 3 Joule/kg °K; 
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c = 1420 m/sec ; a = 2.07 x 10~ 4 / °K; f3 cr = 3 x 10~ 12 / Pascal. 
Inequality (127|) reads 

p«p cr = 2x 10 10 Pascal = 2 x 10 5 atm, 

which is the same critical pressure as in paper [TJ. We observe that Eq. ( 128]) 
yields /3 = 4.98 x 10~ 10 / Pascal which automatically satisfies the inequality 
f l2"B"j) . Liquid water is indeed a very good example of an EQTI liquid. 



4 Quasi-thermal-incompressibility for hyperelastic media 



4-1 Generalities 



As done in paper [TJ and by taking up the definition of pressure presented by 
Flory in 1953 for rubber gum [TBjrB"] and extended for hyperelastic media by 
Gouin and Debieve in 1986 [T7] and Rubin in 1988 [IB] , we can extend the 
results from fluids to thermo-elastic materials. With this aim we define 

C = — -C or C = J^C, (29) 

(detC)s 

where C = F T F is the right Cauchy stress deformation tensor, F is the 
deformation gradient, J = det F = p /p and p is the reference density. The 
specific free energy can be expressed in the form: 

1> = f{p,C,T) , 

where the independent variables p and C are used instead of C. Since det C 
=1, the variable p corresponds to the change of volume, while the tensorial 
variable C is associated with the distortion of the medium: we call C the 
pure deformation of the hyperelastic body. This point is fundamental for the 
decomposition of the stress tensor. If / is independent of C, then we are back 
to the fluid case. It is convenient to introduce the function g such that: 

^ c - r)s/ ( ft (^)i c 4 

Consequently, the free energy of an hyperelastic material can be defined as 
ip = g(p, C, T) where g is a homogeneous function of degree zero with respect 
to C. We deduce the Cauchy stress tensor of the medium in the form |17j : 

t = - p l+ T with p = p 2 ^, T = 2pF^-F T and trr = 0, (30) 

op oC 
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where tr is the trace operator. If g is independent of C, p corresponds to 
the thermodynamic pressure. As proved in |17] . p must be considered as the 
pressure of the hyperelastic material. Let us note that the pressure and the 
change of volume are observable: they can be experimentally measured by 
using a spherical elastic test-apparatus submitted to isotropic stresses. 
The Gibbs equation in the case of elastic materials is [T9|20] : 

TdS = de S -dC , (31) 

where the dot represents the scalar product between matrices and S = 
J F~ l t(F T )~ l is the second Piola-Kirchhoff stress tensor. By inserting Eq. 
(j3DD i into S, Eq. fl3Tl) yields the Gibbs relation: 

TdS = ds + p dV — J~ 3 t • dC, (32) 

with t given by: 

r = — ji F-VF- lT . 
2p 

If we take account of Eq. (129|) . 



and that 



= - p C dV + J§ dC, 
3 



f C = — F^tF-^. C = — tr (f- 1 tF- iT F T f) = — tr r = 0, (33) 



we obtain the Gibbs relation ( 1321) in the final form 

TdS = de + p dV - f ■ dC . (34) 

As in the fluid case, we introduce the chemical potential p = e+p V—TS—t-C. 
Thanks to the orthogonality condition (1331) . p takes the same form as for fluids: 

p = e + pV -TS. 

Equation ( 1341) implies 

dp = Vdp - SdT + t ■ dC. (35) 

Consequently, the change of variables from (V, T) into (p, T) is natural and 
the variable C does not change. Equation fl35|) implies: 



V = p p , S = -p T , t = /ig , e = p - Tp T -pp p , (36) 

with 



dp\ dp\ d V\ 
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Let us assume that the specific volume is written as a function of p,T,C: 

V = V(p,T,C). 
By integration of Eq. (|36|) i . we obtain 

H = f V(p, T, C) dp + fl(T, C), S = -J V T (p, T, C) dp - MT, C) , (37) 

where fx is an additive function depending only on T and C. By substituting 
in Eq. (|36|) /i. we get: 

e = e(T,C) + J V(p, T,C) dp — pV — T J V T dp , (38) 

with the additional function 

e(T,C) = ti-Tth or equivalently fl(T,C) — —T [ ^-^-dT+k{C), (39) 

J T 2 

where k(C) is an arbitrary function of C. The additional function is in the 
same form than for fluids. Therefore 

Statement 4 - For any constitutive functions V = V(p,T,C) , and e 
i : 7 . C), the entropy principle is satisfied if the chemical potential, the entropy 
density and the internal energy are given by Eqs. ( l3~T]) i . f l37|) 9 , f )38|) with Eq. 
(1391)9. respectively. 

In nonlinear theories, the concavity of the entropy density may be valid only 
in some domain of state variables. In particular this is the case in nonlinear 
elasticity where the concavity is in contradiction with the objectivity principle 
if the deformation is large [19] (see also [211122] ) . The concavity of the chemical 
potential is valid only for sufficient small deformations in the neighborhood 
of undeformed configuration. In the decomposition of the stress tensor given 
by Eq. fl30|) . the pressure p is the main observable variable and the pure 
deformation C does not affect the change of variables between internal energy 
and chemical potential. Consequently, taking C constant we conclude that 
inequalities ( ITOlfTTl) are necessary conditions for stability. 

4-2 Quasi-Thermal-Incompressible Solids 

In experiments, the specific volume of incompressible solids changes if the 
temperature varies. In the literature it is usually assumed that J = J{T) 
(or equivalently V = V(T)). Therefore, it is natural to define EQTI solids 
similarly to the definition given for fluids. 
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Definition 2 - A compressible solid satisfying thermodynamic conditions as- 
sociated with the entropy principle and stability is called an Extended- Quasi- 
Thermal-Incompressible solid if there exist V(T) and e(T, C) such that 

V(p, T, C) = V{T)+0{5 2 ) with % = 0(5) and e(p, T, C) = e(T, C)+0{5 2 



Assuming for solids the linear expansion (1241) for V, we derive a similar result 
as for fluids: 

Statement 5 - A thermoelastic material with constitutive equation ([21]) is a 
stable EQTI and tends to an incompressible material under the conditions 

V 1 C P T 

(5 > (5 cr with (3 cr = S 2 , and p <C p cr , with p c 



C P T ' r """" r " 5 k ' 

where V and 5 = «T C 1 are positive constants, while C P may depend on 
the pure deformation C. 



4-3 Case of pure gum rubber 



As in paper [lj, we consider the case of rubber as example of hyperelastic 
material. To verify the conditions EQTI, we get experimental values from the 
literature: numerical values for rubber are obtained in [T^f2"5] : 

p = 10 5 Pascal; C P = 1.9 x 10 3 Joule/kg °K; V = 1.08 x l(T 3 m 3 /kg; 

T = 273 °K; c = 54m/sec; a = 7 x 1(T 3 / °K; = 4.5 x 1(T 7 / Pasca ±; 

(3 cr = 7.6 x 10 -9 / Pascal; p cr = 2.5 x 10 8 Pascal = 2.5 x 10 3 atm. 

We can immediately verify that rubber is also a good example of an EQTI 
medium. Let us note that the considered values for rubber are a little different 
from those in pQ; this is due to the fact that the sound velocity was measured 
at 273°K and not at 325°K, as in [1]. 



5 Conclusion 



Our goal was to present a model expressing the limit case of incompress- 
ible materials and fitting with the principles and conditions of thermodynam- 
ics. Earlier models were thermodynamically deficient. Here, we present a new 
model (which we call Extended- Quasi-Thermal-Incompressible model) start- 
ing from two simple requirements: 

(i) An incompressible medium does not physically exist but it is the limit case 
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of a compressible medium that verifies all thermodynamic conditions. 

(ii) On physical grounds, the volume changes little with the temperature and 

does not practically change with the pressure. 

As result, the conditions for the pressure are the same than in previous pa- 
per [1] but an additional condition associated with the compressibility factor 
must be verified. We point out that a decomposition of the stress tensor for 
hyperelastic media into a pressure term associated with the change of volume 
and a null-trace part associated with the pure deformation of the medium 
leads to analogous conditions of EQTI for fluids and hyperelastic materials. 
Let us note that on a phenomenological basis a theory where the elasticity 
parameters depend on the pressure has been developed in [2~4"f25] . 

From a mathematical standpoint, the difficult question of proving that the 
solutions depend continuously on 5 and that the limit exists as 5 tends to 
zero, remains an open problem. 

Acknowledgments: The authors are grateful to Professor Constantine Dafermos, 
Professor Giuseppe Saccomandi and the anonymous reviewers for their interest on 
this paper and their useful suggestions. 
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